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1. Introduction 

HERWIG is a general-purpose Monte Carlo (MC) event generator for high-energy pro- 
cesses, providing a full simulation of hard lepton-lepton, lepton-hadron and hadron-hadron 
scattering and soft hadron-hadron collisions in a single package, comprising: 

• Initial- and final-state QCD jet evolution with soft gluon interference taken into 
account via angular ordering; 

• Colour coherence of (initial and final) partons in all hard subprocesses, including the 
production and decay of heavy quarks and SUSY particles; 

• Spin correlations in the decay of heavy fermions; 

• Lepton beam polarisation in selected processes; 

• Azimuthal correlations within and between jets due to gluon interference and polar- 
isation; 

• A cluster model for jet hadronisation based on non-perturbative gluon splitting, and 
a similar cluster model for soft and underlying hadronic events; 



A space-time picture of event development, from parton showers to hadronic decays, 
with an optional colour rearrangement model based on space-time structure. 



Several of these features were already present in HERWIG versions 5.1-5.9 and were 
described accordingly in some detail in Ref. [1]. The HERWIG source codes, together with 
other useful files and information, can be obtained from the following web site: 



http://hepwww.ii.ac.uk/theory/seymour/herwig/ 
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2. The MSSM in HERWIG 

Starting from version 6.1 [2], HERWIG includes the production and decay of (s)particles, 
as given by the MSSM. A detailed description of their implementation can be found in [3], 
which complements the current HERWIG manual [4]. The last two successive versions are 
documented in [5], 

The HERWIG particle content is listed in Tab. 1 . For sparticles that mix, the subscripts 
label the mass eigenstates in the ascending order of mass. The two Higgs Doublet Model 
(2HDM) Higgs sector, intrinsic to the MSSM, is also included. The three neutral Higgs 
bosons are denoted by h°, H° and A ( \ whereas the charged ones by H ± . 



Table 1. The MSSM particle spectrum implemented in HERWIG. 



Particle 




Spin Particle 




Spin 


quark 


q 


1/2 


squarks 







charged lepton 


l 


1/2 


charged sleptons 


£l,r 





neutrino 


V 


1/2 


sneutrino 


V 





gluon 


9 


1 


gluino 


I 


1/2 


photon 


7 


1 


photino 


7 


1/2 


neutral gauge boson 


Z° 


1 


zino 


Z 


1/2 


neutral Higgs bosons h°, H°, A 





neutral Higgsinos 


#1,2 


1/2 


charged gauge boson 


W ± 


1 


wino 


W ± 


1/2 


charged Higgs boson 







charged Higgsino 


H ± 


1/2 


graviton 


G 


2 


gravitino 


G 


3/2 



W ± , H mix to form 2 chargino mass eigenstates Xi, xf 
7, Z, Hi 2 mix to form 4 neutralino mass eigenstates X?, X2, 33, X4 
tz,,tjt (and similarly b, r) mix to form the mass eigenstates t\, £2 

HERWIG does not contain any built-in models for Supersymmetry-breaking scenarios. 
In all cases the general MSSM particle spectrum and decay tables must be provided just like 
those for any other particles. A package, ISAWIG, has been created to work with ISAJET 
[6] to produce a file containing the SUSY particle masses, lifetimes, couplings and mixing 
parameters. This package takes the outputs of the ISAIET MSSM programs and produces 
a data file in a format that can be read into HERWIG for the subsequent process generation. 
The user can produce a similar file provided that the correct format is used. To this end, 
we invite the consultation of the ISAWIG webpage: 

tittp://www-thphys. physics. ox. ac.uk/users/PeterRichardson/HERWIG/isawig.htm] 

where some example of input files can be found (the SUSY benchmark points recommended 
in Ref. [7] are already available). Our MSSM conventions are described in [3]. The mass 
spectrum and decay modes, being read from input files, are completely general, however, 
the following caveats are to be borne in mind: (i) SUSY particles do not radiate (which is 
reasonable if their decay lifetimes are much shorter than the QCD confinement scale); (ii) 
CP-violating SUSY phases are not included. 

In addition to the decay modes implemented in the ISAJET package ISAWIG also in- 
cludes the calculation of all 2-body squark/slepton and 3-body gaugino/gluino R-parity 
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violating (RPV) decay modes (alas, RPV lepton-gaugino and slepton-Higgs mixing is not 
considered). 

The emulation of RPV processes is also a feature of the production stage. Tab. 2 
illustrates all MSSM modes available at present (version 6.4). Here, IP ROC is the input 
label selecting the hard process. 



Table 2. The MSSM hard scattering processes implemented in HERWIG. 



IPROC 


2^ 2 MSSM processes in £ + T (I = e, fx) 


700-99 


R-parity conserving SUSY processes 


700 


vr 


— ► 2-sparticle processes (sum of 710— 760) 


710 


t + r 


— ► neutralino pairs (all neutralinos) 


706+4IN1+IN2 


i + r 


— ► XiniXin2 (INI , 2=neutralino mass eigenstate) 


730 


i+r 


— ► chargino pairs (all charginos) 


728+2IC1+IC2 


e + r 


— ► Xici^lC2 ( IC1 ' 2=chargino mass eigenstate) 


740 


£+£- 


— ► slepton pairs (all flavours) 


736+5 1 L 


£+£- 


-> JlM*l,r (IL = 1,2,3 for 1= e, ft, f) 


737+5 1 L 


£+£- 


—> £ L £* L (IL as above) 


738+5 1 L 


£ + r 


— > £l£*r (IL as above) 


739+5 1 L 


£+£- 


— ► £r£* r (IL as above) 


740+5 1 L 


£+£- 


-> vlvI (IL = 1,2,3 for^,^, v T ) 


760 


£+£- 


— > squark pairs (all flavours) 


757+4IQ 


£+£- 


qL,nql,R (IQ = 1...6forg = d...t) 


75 8+4 1 Q 


£+£- 


—> qLq* L (IQ as above) 


759+4IQ 


£+£- 


—> (LQ as above) 


760+4 I Q 


£+£- 


QrQr. (LQ as above) 



Table 2. Continues. 



IPROC 


2^ 2 MSSM processes in £ + T (£ = e, fi) 


800-99 


R-parity violating SUSY processes 


800 


Single sparticle production, sum of 810-840 


810 


£ + £~ — > x°^i, (all neutralinos) 


810+IN 


£ + £~ — > X[ N ^, (IN=neutralino mass state) 


820 


£ + £~ — ► X~ e t (all charginos) 


820+ 1 C 


£ + £~ — > Xic e i~> (IC=chargino mass state) 


830 


— ► and£ + £- -> IfW~ 


840 


£ + T -> u t h°/H°/A Q md£ + r -> £+ H~ 


850 


£+£- -> 


860 


Sum of 870 and 880 


870 


£+r -> £+£~, via LLE only 


867+3IL1+IL2 


£ + r -> ^££2 (ILl, 2 = 1,2,3 for e, ,u, r) 


880 


dd, via LLE and LQD 


877+3IQ1 + IQ2 


£ + r -> d IL id IL2 (IQ1, 2=1,2,3 ford,s,fe) 


Table 2. Continues. 



IPROC 



3000-999 
3000 
3010 
3020 
3030 



2— > 2 MSSM processes in hadron-hadron 



R-parity conserving SUSY processes 

2-parton — * 2-sparticle processes (sum of those below) 

2-parton — > 2-sparton processes 

2-parton — > 2-gaugino processes 

2-parton — > 2-slepton processes 
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Table 2. Continues. 



IP ROC 


2^ 2 MSSM processes in hadron-hadron 


a f\r\r\ r\c\ 

4000-99 


R-panty violating SUSY processes via LQD 


/innn 
4UUU 


single sparticle production, sum of 4010— 4050 




rTi . /-7, v, / A ■ A ^ i. it- n£»ii ("i-qI innc^ 

LljUk ' X, ^% ' ^J^fc X ^ ^"-11 IlCULIallllUa^ 




.7"/ . A , t / A . ^7 -. ^ t / - ^ T T\T — n at i ("rolitir* m o o c ctoto^ 

UjQ>k — ? Xin i ' U'jG'k — ' ain 1 U^-ncULrdiiiiL) mass siaie^ 


4090 


]~i . /-J , i, -\/ j i . A . A , i -\/ yd ('■all r"li Qro'inne^ 

LljUk * X LI j Ilk ' X v"- 11 LllUiglllUS^ 


4020+ I C 


fi'/j, — v „ // ; HacIi — > v ( TP=rharp"inn mass statp^ 

u/j /vjf^ ^ ' J /v A.T{^ i V ^ l cii ^±iikj iiici^ij j lull j 




UjCLk — > Zj ,UjCLk — > l^iVV and CLj(lk — > >C^ VV 


4050 


Ujdk — > £Jh /H /A , Ujd k — > and djd fc — > 


4060 


O J ATA 1 J AHA 

Sum of 4070 and 4080 


4070 


itjfflfe — > u;d m and djdfc — > d;d m , via LQD only 


4080 


itjfflfe — > ^jt fc and djdk — > t^'fc • vla LQD and LLE 


4100-99 


R-parity violating SUSY processes via UDD 


4100 


single sparticle production, sum of 4110-4150 


4110 


~~ *" X "fe> «j'«fc - *■ X w * neutralinos) 


4110 +IN 


Uidj —> Xin"*' — * XiN Ui ( IN as a D° ve ) 


4120 


Uidj — » x u k, djd k —> X di (all cnarginos) 


4120 +IC 


Midj — > Xic" fe ' ^J 6 '*: — * Xic^ ( IC as above) 


4130 


Uidj — > gdk, djdk — > 5«i 


4140 


Ul d 3 fct^ , djdfc t*iZ°, u.dj t*W+ and d 3 d k -» 


4150 


Ul d 3 d* kl h°/H"/A°, djd k u" n h /H"/A", mdj -> djdfc ~* " 


4160 


Uidj — > uid m , djdk — * did m via UDD. 



Table 2. Continues. 



IP ROC 


2 -> 


2 MSSM Higgs processes in £ + £ (£ = e, //) 


910 




" -> j/ e P e /i u + e + e"/i u 


920 


^+^- 


■ Ve» e H" + e+e-H" 


960 




- - 


970 






955 




■ -f 


965 






965 







Table 2. Continues. 


IPROC 


2^2 MSSM Higgs processes in hadron-hadron 


3310,3315 


qq' -> W ± ti',H ± ti> 


3320,3325 


qq' -^W ± H",H ± H" 


3335 


qq' -> H+A" 


3350 


qq -» ^±77^ 


3355 


qg -> H ± H 1L 


3360,3365 


qq^ Z°h ,A°h 


3370,3375 


qq^Z H ,A°H 


3410 


bg — » & ft" + ch. conj. 


3420 


&<? — > 6 + ch. conj. 


3430 


— > 6 ^4° + ch. conj. 


3450 


bg ^ t H~ + ch. conj. 


3610 




3620 


qq/gg -> -ff 


3630 
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Table 2. Continues. (For the definition of ISQ, see [5].) 



IPROC 


2 — > 3 MSSM Higgs processes in hadron-hadron 




99 /qq — » qq 


3200+1 SQ 


gg/qq -> gg *h,H,A 


3500 


bq — * bq H ± + ch. conj. 


3710 


/ — 7T7T 

gg -» g'g' /i u 


3720 


gg -» g'g'tf 


3810+IQ 


gg + qq ^ QQh u (IQ for Q flavour) 


3820+IQ 


33 + gg -> QQH° (") 


3830+IQ 


53 + gg -> QQA (") 


3839 


<7g + gg — > + ch. conjg. 


3840+IQ 


gg — * QQh° (IQ as above) 


3850+IQ 


gg -» QQ//° (") 


3860+IQ 


-» QQA° (") 


3869 


gg — » btH + + ch. conjg. 


3 870+ I Q 


gg — > QQh° (IQ as above) 


3880+IQ 


gg QQi? (") 


3890+IQ 


gg ^ Q_QA° (") 


3899 


gg — > btH + + ch. conjg. 



The Matrix Elements (MEs) used for the hard scatterings have been listed in [3], with 
the treatment of the different colour connections in SUSY QCD processes described in [8] 
(which also superseeds the standard QCD algorithm of [1]). 

Spin correlations are available in processes where SUSY particles are produced (and top 
quarks and r leptons as well), as described in [9]. Whenever these particles decay, 3- and 
4-body MEs are used to describe their dynamics (with or without the spin correlations). In 
the case of r leptons, an interface to TAUOLA [10] is also available. 

Finally, polarisation has been implemented for incoming leptonic beams in SUSY pro- 
cesses. These effects are included both in the production of SUSY particles and via the spin 
correlation algorithm in their decays. 
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